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Abstract 

We present fully dynamic simulations of heavy ion collisions at RHIC energies within the perturbative QCD-based 
partonic transport model BAMPS, focusing on the simultaneous investigation of jet-quenching and elliptic flow. The 
model consistently features elastic and inelastic 2 <-> 3 processes, the latter being based on the Gunion-Bertsch matrix 
element. We discuss first attempts to extend the model to include light quark degrees of freedom and study the energy 
loss of high energy gluons and quarks in a static partonic medium. The difference between gluons and quarks in 
inelastic processes is found to be weaker than expected from color factors, due to a self-quenching effect associated 
with a cut-off modeling the LPM effect. 
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1. Introduction 

Two of the most important phenomena observed in collisions of heavy nuclei at the Relativistic Heavy Ion Collider 
(RHIC) go by the names of jet quenching and elliptic flow. 

The RHIC experiments have established that particles with high transverse momenta are suppressed in heavy ion 
collisions with respect to a scaled proton-proton reference (HE). This quenching of jets is commonly attributed to 
energy loss on the partonic level. It could therefore provide means of investigating the properties of the medium, 
the quark-gluon plasma (QGP), that causes the modifications of the high-p T particles as they traverse it after their 
production in initial hard processes. Jet quenching on the single hadron spectra level is usually quantified in terms of 
the nuclear modification factor Raa- 

The collective flow of the created medium 0) is usually quantified in terms of the coefficient V2 in a Fourier 
expansion of the angular distribution dN/d(p and in this context often referred to as elliptic flow. The comparison 
of hydrodynamic calculations to data indicates that the viscosity of the QGP is quite small Q, possibly close to the 
conjectured lower bound n/s = l/(47r) from an AdS/CFT correspondence (6). 

The energy loss of partonic jets on the other hand can be treated in terms of perturbative QCD (pQCD) and most 
theoretical schemes attribute the main contribution to partonic energy loss to radiative processes ITl-HOll. 

It is a major challenge to combine high-p T physics and bulk evolution within a common framework. Recently, 
efforts combining pQCD-based energy loss calculations with hydrodynamic modeling of the medium have been in- 
tensified, using results from hydrodynamical simulations as an input for the medium evolution in jet-quenching cal- 
culations (see ifTTTl for an overview). However, these approaches still treat medium physics and jet physics in the QGP 
on different grounds. No schemes are available that cover the full dynamics of the interplay between jets and the 
medium, i.e. that consistently include modifications of the medium caused by the traversing jets. 

Partonic transport models might provide means to investigate bulk properties of the QGP and high-energy parton 
jets within a common physical framework automatically including the full dynamics of the evolution of the system. In 
previous publications E2l[T3l we have explored the capabilities of the transport model BAMPS (Boltzmann Approach 
to Multi-Parton Scatterings) with this goal in mind. 
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2. The transport model 



The microscopic transport model B AMPS (HI Q21 is a tool to simulate the evolution of the QGP stage of heavy 
ion collisions. Partons are considered as massless Boltzmann particles whose interactions are based on leading order 
pQCD matrix elements and consistently include creation (2 — > 3) and annihilation (3 — » 2) processes. The test particle 
method is introduced to reduce statistical fluctuations. 

BAMPS has so far been restricted to gluonic degrees of freedom, that already allows for the investigation of 
crucial features of the medium evolution in heavy ion collisions. Nevertheless, the contribution of light quarks is of 
great interest and we will discuss a possible extension of BAMPS to quark degrees of freedoms. 

For elastic interactions Debye screened cross section in small angle approximation are used, for example 

dcr gg ^ gg 9na 2 dcr qg ^ qg 2na 2 



dq\ (q 2 ± + m 2 D ) 2 dq\ (q 2 ± + m 2 D ) 2 

The Debye screening mass is computed via m 2 D = dc7ia s J ^3];(N c f g + Nff q ) from the local gluon and quark 
distributions, f g and f q , where do = 16 is the gluon degeneracy factor for N c = 3. 




Figure 1: Left panel: Total dE/dx for a high energy gluon and a high energy light quark as a function of energy. Medium temperature T = 0.4 GeV. 
Right panel: Ratio of the mean free path for gluons to the mean free path for light quarks for purely 2 — > 2 interactions and for 2 — > 3 processes. 

Inelastic processes are treated via the Gunion-Bertsch matrix element lfT6l . For gg —> ggg it reads: 
■ ,2 _ 12n 2 a 2 s 2 4Sna s q 2 ± / v 

|At ™! = (qi + m^ k 2 ± [(^- q J 2 +m 2 D ] ^ " ' (2) 

where q ± and k ± denote the perpendicular components of the momentum transfer and of the radiated gluon mo- 
mentum in the center of momentum (CM) frame of the colliding particles, respectively. The principle of detailed 
balance yields the matrix element for the back reaction. Inelastic processes involving quarks, e.g. qq — > qqg, are 
dealt with by considering that the 2 — > 3 matrix element factorizes into a probability for the emission of a gluon 
and an elastic contribution. This in mind, the matrix element for e.g. qq — > qqg can be obtained from ^ via 

\M qq ^ qqg \ 2 = \ Mqq ^ qq \ \M gg ^ ggg \ 2 . Taking the small angle approximation for the elastic matrix elements, the prefactor 

i i 2 i i 2 

\M qq ^ qq \ I \M gg ^g g \ simply reduces to a momentum independent color factor. 

The Theta function in ^ is an effective implementation of the LPM (Landau, Pomeranchuk, Migdal) effect ifTTl . 
that describes coherence effects in multiple bremsstrahlung processes. This interference effect cannot be incorporated 
directly into a semi-classical microscopic transport model such as BAMPS, so the cut-off (A - r) ensures that 
successive 2 — > 3 processes are independent of each other, r is the formation time of the gluon emitted with transverse 
momentum k ± and A denotes the mean free path, i.e. the time between successive interactions, of the parent parton. 
When comparing the formation time to the mean free path of the parent parton, attention needs to be paid to the 
reference frames. This renders the cut-off dependent on the boost J3 between the plasma rest frame and the center of 
momentum (CM) frame in which ^ is evaluated fT3ll and numerically further complicates the calculations. 



3. Results 



Studying the evolution of high energy partons in a static thermal medium with fixed temperature, we find that 
within our setup the mean energy loss per unit path length dE/dx is strongly dominated by inelastic 2 — > 3 processes. 
dE/dx in a static medium is calculated from the average energy loss per collision and the average interaction rate as 
dE/dx = dE/d(ct) = Yji(AE l )R\ where the sum runs over all possible interaction processes. 

While dE/dx from binary processes follows the well know T 2 ln(ET) dependence and is moderate at dE/dx gg ^ gg « 
1.2GeV/fm for E = 50GeV and T = 0.4 GeV, the energy loss from interactions according to the Gunion-Bertsch 
matrix element rises roughly linearly with E and is rather strong with dE/dx gg ^ ggg « 32.6GeV/fm for the same 
medium parameters fT3ll . The total dE/dx can be seen in the left-hand panel of fig. [T] 

For elastic processes the mean free path and dE/dx of a high energy gluon and a light quark differ by roughly a 
factor of 2 as expected from the color factors. For 2 — > 3 processes however, the difference is much weaker. Only 
about 20% in the energy loss dE/dx and the mean free paths, see fig. [T] This is due to the cut-off (A - r) that 
basically self-quenches the difference in the matrix elements for gluons and quarks. When the mean free path is 
larger, as is expected for quarks from the color factors, the step function for the integration over the matrix element 
^ becomes less restrictive, yielding a larger cross section and thus reducing the mean free path. The results shown 
in fig. [I] are self-consistently obtained from iterative calculations of the mean free path A. 
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Figure 2: Left panel: Elliptic flow V2 for gluons in simulated Au + Au collisions at 200 AGeV with b = 7 fm. 

Right panel: Gluonic Raa as extracted from BAMPS simulations for b = fm and b = 7 fm. For comparison experimental results from PHENIX 
1 18 1 for tt° are shown for central (0% - 10%) and off-central (20% - 30%) collisions. See 1 13 1 for details. 

Since the difference in the differential energy loss in static media from our treatment of quark inelastic interactions 
is rather weak, only a small difference in the nuclear modification factors of gluons and light quarks is to be expected 
in fully dynamic simulations of Au + Au collisions at RHIC energies. However, these calculations have not been fully 
completed yet and will be presented in a later publication. 

For a purely gluonic medium Raa in central Au + Au collisions (initial gluon distribution from a Glauber plus 
mini-jet model) is flat at roughly ^f^° ns ~ 0.053 and in reasonable agreement with analytic results from Wicks et al. 
(l0l[T2]|, though the suppression of gluon jets in BAMPS appears to be slightly stronger. In non-central collisions the 
suppression is reduced but Raa remains flat. The gluon Raa at b = 7 fm increases by a factor of roughly 2.5 compared 
to b = Ofm. The relative increase is less pronounced (~ 1.7) in the pion data, comparing 0%-10% and 20%-30% 
central collisions. See right panel of fig. [2] The V2(pr) of high-p T gluons at b = 7fm is rising up to V2 ~ 0.12 at 
p T ~ 4 GeV and slightly decreases afterwards (left panel of fig. [2]). This behavior is in good qualitative agreement 
with recent RHIC data fT9l 



4. Summary 

The transport model BAMPS can be used to study the energy loss of high-p T partons and the elliptic flow within 
a common framework. We find gluonic Raa for central Au + Au to be in reasonable agreement with analytic results 
and the integrated elliptic flow can be reproduced nicely [ 20 1. Inelastic processes based on the Gunion-Bertsch matrix 
element strongly dominate the energy loss of high-p T partons. 
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First attempts to extend the model to include light quarks show that the energy loss of quarks and gluons differs by 
less than the expected color factor due to a self-quenching effect caused by the effective implementation of the LPM 
suppression. This will be studied in more detail in upcoming works, together with the application of fragmentation 
functions to the high-p T sector allowing for more direct comparison to hadronic observables. Also the simulation 
of LHC energies and a careful scan of various impact parameters is on the agenda, improving the comparability to 
experimental data. 
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